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The static color-Coulomb interaction potential is calculated as the solution of a non-linear inte- 
gral equation which arises in the Hamiltonian Coulomb gauge when the restriction to the interior of 
the Gribov horizon is implemented. The potential obtained is in qualitative agreement with expec- 
tations, being Coulombic with logarithmic corrections at short range, and confining at long range. 
The values obtained for the string tension and Njjg are in semi-quantitative agreement with lattice 
Monte Carlo and phenomenological determinations. 
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t-H ' The interaction energy E(r) of a pair of heavy quarks at separation r is a prominent feature of QCD phenomenology. 
;_i From fits to the charmonium and bottonium spectra E(r) is well-known in the range 1 (GeV) -1 < r < 4 (GeV) -1 . 
, Wj ' To calculate E(r) from first principles, one makes a Born-Oppenheimer approximation, in which heavy quarks move 
, slowly, while gluons and light quarks rapidly adjust to the instantaneous position of the heavy quarks. In this 
approximation, E(r) is the potential energy that appears in a non-relativistic Schroedinger equation that describes 
| the slow motion of the heavy quarks, and from which the phenomenological fit for E(r) is obtained. It is calculated 
as the gauge-invariant ground-state energy = E(r)^Q, where Ti is the field-theoretic QCD Hamiltonian that 

describes the dynamics of gluons and light quarks in the presence of an external quark-anti-quark pair at fixed 
. separation r. In the absence of light or dynamical quarks, it is generally believed that E(r) grows linearly at large r, 
Tj" ' E(r) ~ Kr, with string-tension if, and this belief is supported by lattice-gauge calculations Pair production of 
C^l , dynamical quarks from the vacuum causes "breaking" of the string at large r. However, in the above energy range 
■ where E{r) is known phenomenologically, string-breaking is not yet manifest, and we shall ignore dynamical quarks 
in the following. 

The Coulomb-gauge is a so-called "physical" gauge in which the conjugate dynamical variables are the 3- 
[ dimensionally transverse fields A — A tr and (—E tr ), where E^ { = (-E* r )S.i ~ (Viy )^ and <p is the color-Coulomb 
operator. Although the Coulomb-gauge is not known to be renormalizable, a lattice Coulomb-gauge Hamiltonian has 
been derived recently from the transfer matrix of Wilson's gauge-invariant Euclidean lattice gauge theory so that 
gauge-invariant quantities such as E(r) have the same continuum limit as in Wilson's theory. For brevity, equations 
will be written formally in continuum notation. The exact lattice equations may be found in Q . 

The transversality condition on A is not a complete gauge fixing because of the existence of Gribov copies ||. 
To obtain a complete gauge fixing, we adopt the minimal Coulomb gauge. In this gauge, the set A of gauge-fixed 
configurations, called the fundamental modular region (FMR), consists of configurations A for which the Hilbert norm 
\\A\\ 2 = J d 3 x\A\ 2 of A is minimized with respect to local gauge transformations g: A = {A : \\A\\ < \\A 9 \\ forallg}. 
Here A 9 = g^Ag + g^Vg is the gauge transform of A. Configurations with this property are transverse, A = A tr — so 
this gauge falls into the class of Coulomb-gauges — and moreover the Faddeev-Popov (FP) operator M(A) = —D ■ V 
is symmetric M(A) = M^(A) and positive M{A) > for all A S A [Q. Here D = D(A) is the gauge-covariant 
derivative D a > c (A) = 5 a ' c V + g f abc A b . In the following we write A for A tr . 

With E = E tr — V(f>, the non-Abelian Gauss's law constraint, D ■ E — gop qu , may be written M<p = g$p. Here p qu 
is the color-charge density of the quarks, and p a = —f abc A tr - b E tr ' c + p a qu is the color-charge density of the dynamical 
degrees of freedom. Because M(A) is a positive operator in the minimal Coulomb gauge, the color-Coulomb field 
<j) may be expressed unambiguously in terms of the dynamical degrees of freedom, <f> — 5o^ _1 P, and so also the 
color-electric field E — E tr — V(f>. For wave functionals 3>(j4) and ^(A) defined for A in A, and with E tr — iS/8A tr , 
the Coulomb-gauge Hamiltonian is defined by the quadratic form 

($,ii cou i*) = J dA tr a J d 3 x2- 1 [(E?<P)*E°-$> + <i>*B?B^} , (1) 
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where E has just been defined, Bf = VjA k - V k A 3 + f abc A b A c k for and k cyclic, a = det[M(A)/M(0)], and 
the quark Hamiltonian is suppressed. This Hamiltonian is symmetric with respect to the inner product ($, 4") = 
J A dA tr a^*^, and a(A) is positive because M(A) is. 

The main novel ingredient in the present approach is the implementation of the restriction of the preceding integrals 
to A, a, non-perturbative effect [Q. For finite quantization- volume V, the exact boundary of A is not known. However, 
in y it is argued that, for periodic boundary conditions and in the limit of large V, this region is adequately 
approximated by G/V < 0, and moreover that this restriction may be implemented by use of the effective Hamiltonian 
H CO ui =^ H c ff = H cou i + 70G. Here the horizon function G(A) is defined for structure group SU(N) by 

G(A) = f d 3 x d 3 y D a £D a J{M- l y> d {x, y) - 3(N 2 - l)V . (2) 



The term jqG makes the wave-functionals in the Fock space of H c g vanish rapidly outside A, so the restriction on the 
integrals may be ignored, J. dA tr => J dA tr . Here 70 is a thermodynamic parameter that sets the scale for hadronic 
masses. Its value is determined by the horizon condition (G) /V — 0. Here and below, the expectation- value is taken in 
the ground state of H e g . The horizon condition expresses the fact that, for large V, the probability gets concentrated 
on the boundary G/V = 0. 
The Coulomb energy 

Seoul = [dA tr a [ d 3 x2- 1 \g V(M- 1 p)<f>\ 2 (3) 



contributes additively to the energy of a state <&, the remaining terms being positive. Consider the contribution to 
£"coui from the part of M _1 which is diagonal in momentum space, 

V-\M- l ) q ^ h = fd 3 ye lq y [V- 1 [d :i .r:M ; i,.„ : ,. ;/ ,/ . (4) 



For each y, the term in brackets is a bulk quantity per unit volume. By translation invariance, its covariance matrix 
is of order V^ 1 . Consequently, for large V, it approaches its mean-field value 

r^M- 1 ),,^ = c( q )s a - b + oiv- 1 ' 2 ) , (5) 

where 

C(q)5 a * = [ d 3 ye i ^((M- 1 ) x , a;x+a , 6 ) (6) 



is the 3-dimensional FP propagator. Thus <p and Seoul receive the contributions 4>{q) = 9oG(q)p(q) and 

J dA tr a (2V)- 1 E g q 2 \goC(q)p(qm 2 , (7) 

where the color-charge density operator p(q) satisfies \m\ q ^p a (q) = Q a . Here Q a , with [Q a ,Q b ] = f abc Q c , is the 
total color-charge. 

The FP propagator is of the form C~ l (q) = q 2 — qiTiijqj, where Ejj(g) is a self energy. The horizon condition is 
equivalent Q to the condition Eij(0) = This gives C- 1 (q) = gj[S,-j(0) — E^g)]^-, so C(q) is more singular at 
<7 = than q~ 2 . We conclude that the restriction to the FMR produces a long-range color-Coulomb field, as originally 
foreseen by Gribov || . 

The self-energy E satisfies a Schwinger-Dyson equation which we write symbolically qT,q = R(q) = g 2 . J T D C T, 
which gives C^ 1 (q) = R(Q) — R(q). Here D is the gluon propagator, and T and To are the exact and zeroth-order 
ghost-ghost-gluon vertex functions. When R is evaluated in the ground state of H e g, this equation || provides the 
mean- field self-consistency condition (MSC) that determines the mean- field function C(q). 

So far our results are exact, but in the following we resort to approximations to solve the MSC. This equation 
contradicts the usual perturbative expansion, as one sees from the appearance of (go) 2 in R, whereas C, D, and T 
are nominally of leading order (go) - To obtain an expansion in powers of go that is consistent with the MSC, we 
assume that C(q) is of leading order (go) 1 , namely C(q) = g _1 u(g) + 0(1), and that all other correlation functions 
are analytic in g . To leading order, the MSC is an integral equation that determines u(q): 

u-\q) = iV(27r)- 3 / d 3 kDo(k)[q 2 -(q-k) 2 /k 2 ][u(k)-u(k-q)] . (8) 
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The kernel Dq is the gluon propagator 

D (k)(6 itj fc^/fc 2 )^ = J d 3 X e- ik -*{A b x+y .A c yd ) , (9) 
evaluated in the ground state of 

H Q = J2 ^ «l,A«fe,A + (Wy 1 P-kVkP a k , (10) 

fe,A k 

where a\ x and a^^x are creation and annihilation operators for A tr and E tr , A is a two- valued polarization index, and 

ti>k = (fc 2 + /i 4 ?/^) 1 / 2 . This is the zeroth-order part of H c g , obtained by systematically expanding H e g in powers of c/Oj 
with 70 scaled according to 70 = (27V go) 1 1 1 ^ ■ I n addition to a harmonic oscillator Hamiltonian, Hq also contains a 
color-Coulomb interaction Hamiltonian with interaction potential Vk = fc 2 u^ that is independent of 30, as is consistent 
with dimensional transmutation. 

This term prevents us from calculating Do(k) exactly. We neglect it, so Dq is approximated by D^ °\k) = (2u>k)^ 1 , 
and we call w q = u q ^ the solution to (||) with kernel . We expect that Do(k) is even more suppressed at low k 

than D { °\k), so by (|) u q is enhanced at low q compared to w q corresponding to an even longer range force than we 
find. 

The asymptotic form of w q at high and low q has been determined analytically B , and it has been found numerically 
for intermediate q with 2% accuracy (see [| 10 for details). We write w(q) = g c (q)/q 2 , where g c {q) is a renormalized 



running Coulomb-coupling constant, and we express our results in terms of the static color-Coulomb potential v q — 
q 2 w 2 (q) = g 2 (q)/q 2 , which appears in Hq, and the Fourier transform of v q (normalized for a pair of external quarks), 
which is given (apart from an additive constant) by 

poo 

V(r) = —(N 2 - 1X47V7T 2 )- 1 / d qg 2 c (q) (qr)- 1 sin(qr) , (11) 



corresponding to a force f(r) = —V'(r). 

The asymptotic form of g c at high q is given by 
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= b c t + 6 _1 6 c ln(6 c t) + Oir 1 \nt) , (12) 



where t = \n(q/m), b c = (3tt 2 )~ 1 N, m = Cfi, and C is a constant. The first and second terms are of the form of 
the one- and two-loop contributions to the running coupling constant g Tg of the perturbative renormalization group, 
but the coefficients are different. This difference arises because at high momentum the interaction is not purely static 
Coulombic, and it is verified in Q that the difference in the first coefficient is correctly accounted for by terms that 
are neglected here. The coefficient of the second term is a new result. The limiting behavior of the force /(r) at small 
r is given by 

f(r) w37r(7V 2 -l)(4^ 2 r 2 )- 1 [log(A| s r 2 ) + 3- 1 loglog(A 2 r 2 )]- 1 , (13) 

where = e 7_1 m, 7 is Euler's constant, and A is a constant. 

The confinement properties of the theory in the present approximation are determined by the asymptotic form 
of g c (q) at low q, which is given by g c {q) = B^i/q) 4 / 3 , where B~ 3 / 2 = 7V7r- 2 r(8/3)r(2/3)/r(16/3) (see f§). This 
corresponds to a color-Coulomb potential V(r) ~ r 5 / 3 which rises more rapidly with r than a string tension Kr. 
This somewhat surprising result may be an artifact of the approximations made. On the other hand V{r), which 
appears in the (approximate) quantum-field theoretic Coulomb-gauge Hamiltonian, must be distinguished from the 
gauge-independent quark-pair energy E(r) discussed in the introduction, and it is shown in [Q that E(r) < V(r). 
If V(r) does grow more rapidly than Kr at large r, then the wave-function 1 ^(A tr ) 1 whose defining property is to 
minimize E(r), adjusts itself so E{r) rises no more rapidly than Kr [ pd| . It does so by changing the two superposed 
spherically symmetric long range color-Coulomb fields into a flux tube. 

We shall compare our results for V(r) with phenomenological fits to E(r), to see if there is a range of "small" r for 
which V(r) agrees with E(r), as suggested by asymptotic freedom. To this end, we consider two phenomenological 
models: the Cornell potential [j] and the Richardson potential 0, both of which fit the cc and bb spectra well. 

To connect dimensionless quantities and the real world, we fix the length scale by using Sommer's |l2] | dimensionless 
phenomenological relation r\ f(ro) = —1.65 which holds for the Cornell force at ro a = Rq = 2.48 (GeV) -1 . From 
the value of a in (GeV) -1 we obtain the string tension \/~!r = y/mm r [— f(r) } a -1 = 518 MeV. [By this definition, 
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we are evaluating a where / (r) = 0, i.e where the potential is approximately linear.] It is not easy to estimate an 
uncertainty for the string tension. However, its value seems to depend very weakly on the values of the parameters 
of our trial solution ||. If we identify the parameter = e 7 ~ 1 m in ( p"3| ) with the corresponding physical parameter 
[[l|, then from the relation [||| A-^ = A R exp (-7 + 1 - 31/66) we obtain A^ = 124 ± 12MeV (see ||]§). In 
Figure [l]a we plot our result for f(r), the two Cornell forces |Ej, and the Richardson force B. Our force gets its 
maximum value at a separation of about 2.5(GeV) _1 and it is almost constant up to 4(GeV) , the variation being 
of order 12%. 

These results depend on the phenomenological conditions that we have used to set the length scale. If we use 
the dimensionless relation r% f(r ) — —1.35, which holds for the Richardson force at i?o = 2.48(GeV)~ 1 , we obtain 
\fa = 468 MeV, A-r^j = 118 ± 12MeV and the plot shown in Figure [j]b. In this case the agreement is even better: 
our force reaches its maximum value at a separation of about 2.75(GeV)~ 1 , and its variation at 4(GeV)~ 1 is of order 
8%. 

As an exact result, we have found that the restriction to the fundamental modular region causes an infrared 
singularity of the color-Coulomb propagator and thus a long-range color-Coulomb potential. Although our calculation 
of this potential required possibly severe approximations that are describe above, nevertheless it gives results which 
are in qualitative agreement with phenomenologically determined potentials. The fit is surprisingly good. After 
setting the scale by Sommer's method, the force is in qualitative agreement with phenomenological models, and the 
values obtained for the string tension and Ajjg are in semi-quantitative agreement with lattice Monte Carlo and 
phenomenological determinations (see jj],U and jl5|]). It would appear that the approximate equality V(r) w E(r) 
extends to the range r < 4(GeV) _1 , and moreover that the approximations made in our calculation of V(r) do not 
qualitatively destroy this agreement. Although there is no a priori reason to expect that vacuum polarization of gluons 
should not be important in this range, this may not be so surprising after all. For although vacuum polarization (pair 
production) of quarks does "break" the string, this is not yet manifest for r < 4(GeV)~ 1 . 

We are grateful to Tony Duncan, Martin Schaden and Alberto Sirlin for informative discussions. 
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FIG. 1. Plot of: (i) our force f(R) (the curve which is increasing negative at large R), (ii) the forces (for the cc and the bb 
cases) derived from the Cornell potential (the two curves very close to each other) , and (iii) the Richardson force (the highest 
curve). In case 1 we set the length scale by r$ f(ro) = —1.65 from the Cornell force, while in case 2 we use the condition 
r o f( r o) = —1-35 from the Richardson force, both at R — 2.48(GeV)~ 1 . 
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